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Abstract 


Part of the order a s correction to the inclusive B —* X s 7 photon spectrum is 
enhanced by /rnf). We discuss its origin and sum the corrections pro¬ 

portional to [a s log(ml/fi 2 )] n to all orders. These are the calculable leading 
logarithms in the parton fragmentation functions of a quark or gluon into a 
photon. Although the gluon fragmentation into a photon starts only at order 
a 2 , its contribution is of the same order as the s quark’s in the leading log 
sum. For not too small values of the photon energy, the resummation yields 
a moderate suppression. In the standard model, the coefficient of the opera¬ 
tor whose matrix element gives rise to such terms is small. A measurement 
of the photon spectrum around 1 GeV would provide a theoretically clean 
determination of C%, the Wilson coefficient of the b —> s g operator. 


*Work supported in part by the U.S. Dept, of Energy under Grant no. DE-FG03-92-ER 40701. 
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I. INTRODUCTION 


Inclusive B —> X s 7 decay is very sensitive to physics beyond the standard model (SM) 
|L] ||. As any flavor changing neutral current process, it can only arise at one-loop level 
in the SM, and therefore possible new physics can yield comparable contributions. As the 
recent CLEO measurement [7j] excludes large deviations from the SM, a next-to-leading 
order calculation is essential. The absence of such a calculation will soon become the major 
obstacle in improving the bounds on new physics from this process. Certain parts of this 
computation have already been carried out ||||, but the technically most challenging parts 
are yet to be calculated. 

It was observed recently || that moments of the photon spectrum can be predicted 
to order a s accuracy by knowing the coefficients of the effective Hamiltonian only to the 
presently available leading logarithmic accuracy. This will provide a model independent 
determination of the b quark pole mass (see also [1(J), that is, the A and Ai matrix elements 
of the heavy quark effective theory. 

In this paper we point out that part of the next-to-leading order correction to the photon 
spectrum induced by the b —> sg operator is enhanced by log(m%/m%). We describe these 
contributions in section II. In section III we resum the series of leading logarithms of the form 
[a s log(m^//.i 2 )] ri . We discuss the meaning of the infrared scale /j and how nonperturbative 
phenomena enter the results. The resummation of leading logarithms does not eliminate the 
logarithmic dependence on the b quark mass, which formally cancels the a s suppression of 
this term. In section IV we discuss possible phenomenological implications. 

II. log(m^/m^) TERMS IN THE ORDER a em a s RESULT 

I 11 the standard model, B —> X s 7 decay is mediated by penguin diagrams. QCD cor¬ 
rections to this process form a power series in the parameter a s log(M^/m^), which is too 
large to provide a reliable expansion. Therefore, it is convenient to integrate out the virtual 
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top quark and W boson effects (and possible new physics) at the W scale, and sum up the 
large logarithms using the operator product expansion and the renormalization group. We 
work with the operator basis and effective Hamiltonian of Ref. [2J 
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For the present paper only two of these operators are directly relevant 
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where a, (3 are color indices, and P Ry L = 1(1 ± 75 ). Throughout our discussion C 7 and Cg 
will refer to the “effective” Wilson coefficients [Q that include the leading order contribution 
of the four quark operators. 

We are interested in the order a s correction proportional to Cg, which is the order a em 
corrections to B —> Ah. g. The hnite part of the spectrum is given by 
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Here C F = 4/3 in SU( 3), 
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(2.4) 


is the leading order contribution to the B X s 7 decay rate given by the matrix element 
of Oj, and we introduced the dimensionless parameters 
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x = 


(1 — r) m b ’ 


r = 


9 

mt 


(2.5) 


The variable x corresponds to in the free quark decay model. The extra factor of 

9 suppression in eq. (|2.3|) results from the square of the b or s quark’s electric charge. 

This contribution to the total B —> X s 7 decay rate is infinite. The divergence is canceled 
by the virtual correction to the B —> X s g decay mode (with no photon in the final state). 
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While soft photons could be summed to all orders in a em to make the E 1 —> 0 limit smooth, 
in finite orders of perturbation theory, it only makes sense to talk about the total inclusive 
B X s 7 decay rate with an explicit lower cut-off on the photon energy. 

The more interesting feature of the spectrum in eq. (|2.3|) is that it is not finite in the 
mb/m s —> oo (r —> 0) limit. In particular, the photon spectrum calculated in dimensional 
regularization with vanishing light quark masses would contain a part proportional to 1 /{D — 
4) for any value of E 1 . This is a consequence of the fact that a massless quark cannot be 
distinguished from the same massless quark and one (or more) collincar massless partons 
with the same total charge and energy. The singularity is eliminated once a nonzero quark 
mass is included, but it results in the above logarithmic dependence of the spectrum on r. 

Pictorially, this is simplest to understand by looking at the Dalitz plot in the r —> 0 limit 
(in the f ? 7 — E g variables), shown in Fig. 1. The shaded area represents the available phase 
space, the thick line denotes where collinear singularities occur, and the black box indicates 
where the virtual corrections contribute. Fig. 1/a corresponds to the square of the invariant 
matrix element of 0 7 . In this case both the collinear singularity and the virtual correction 
occur at E 1 = mb/2. For this contribution, there is no problem whatsoever carrying out the 
calculation in the rn s = 0 limit. The reason is that the kinematic variable (1 — x) regularizes 
all singularities, and one finds the usual picture of Sudakov exponentiation near ff 7 = mb/2. 
However, the situation is more complicated for the contribution of the square of the invariant 
matrix element of 0 8 , shown in Fig. 1/b. In this case, collinear singularities occur at any 
value of f? 7 ; these are only canceled by the singular part of the virtual correction in the total 
decay rate, but not in the photon spectrum. 

If the hadronic final state X s formed a jet, by excluding photons within a certain angu¬ 
lar cut around the jet axis, the above discussed collinear singularity would be eliminated. 
However, at CLEO B meson decay products are distributed in angle. Thus, when the ex¬ 
perimental analysis excludes events in which a photon and a charged particle overlap in 
the detector, this cut is not related to the direction of the decay product s quark at short 
distances. 
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III. RESUMMATION 


To identify the origin of the leading logarithms of the form [a s log(m^//z 2 )] n , and to 
resum them to all orders, we observe that the analytic structure and r-dependence of the 


O8 contribution is very similar to the photon structure function, as discussed in [11]. In the 
latter case one is interested in a hard photon scattering off another (almost) on-shell photon, 
while in the present problem an on-shell photon is created in the decay of a heavy particle. 

We work to fixed (first) order in a em and to all orders in a 8 . It is simplest to consider 
the discontinuity of the forward scattering amplitude corresponding to intermediate states 
with an s quark, a photon, and any number of light quarks (u, d, and s ) and gluons. In 
an axial gauge the leading logarithms come from those diagrams in which at order a™ there 
are precisely n ways of separating the external photon legs from the other external legs by 
cutting two internal lines (see, e.g., [0). This implies that the result is given by the sum 
of two sets of ladder diagrams. One corresponds to the decay function of the s quark into a 
photon, and the other to that of the gluon. This second set of ladder diagrams starts only 
at order a 2 sl but also yields leading logs. From these arguments it also follows that only the 
square of the Og operator’s invariant matrix element can yield these type of leading logs. 

Decay functions are typically not calculable from first principles; moreover, they are 
usually divergent in perturbation theory. As the photon only couples to massive partons 
(quarks), in a perturbative calculation the light quark masses provide an infrared cut-off. 
That is why eq. fl2.3|) in the previous section contains a logarithmic dependence on m s . 
Higher order corrections would also include terms proportional to log rn u and logm^. How¬ 
ever, in real QCD there is no such dependence on the light quark masses. Non-perturbative 
decay functions are well-defined, measurable physical quantities. 


The crucial observation of Ref. [lT[] was that for a parton interacting with an on-shell 
photon, the absolute normalization of the leading logarithms can be calculated in pertur¬ 
bative QCD. (In other cases only the Q 2 evolution is calculable.) This large Q 2 (or large 
ml) behavior is determined by the matrix element of the photon operators between the 
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external photon states. The contributions of quark and gluon operators are logarithmically 
suppressed. To sum the leading log ladders, one introduces an infrared scale /x, which is 
appropriately chosen to be a few times Aqcd- The question of the //-dependence is then 
non-leading in logm&. Addressing it consistently involves a combination of perturbative and 
non-perturbative (he., decay functions deduced from other experiments) terms of compara¬ 
ble magnitudes. All strongly coupled long distance effects are absorbed into the definitions 
of the decay functions. 

The resummation of these leading logarithms is important, since mj // implies 
log(m 2 //x 2 ) 1 /a s (rrib). Only in this limit can inclusive decay rates be calculated model- 

independent ly, based on an operator product expansion [Hj. The result of the resummation 
of the leading logarithms is given by the sum of the decay functions of the s quark and the 
gluon into a photon 
dT« 


da; 


resummed 


= T(b -»■ sg) x [D s ^(x) + D 9 ^{x)}. 


(3.1) 


Here T(6 — > sg) = T 0 (a s Cp C\fa em C 2 ) is the tree-level b — > sg decay rate. The calculable 
leading logarithmic part of the decay functions are simplest to write down in terms of their 
moments (n > 2 ) [|H 
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Here V n = f ( ] dxx n ~ 2 [1 + (1 — a;) 2 ], e s = —| is the s quark’s electric charge, and (e 2 ) denotes 

the average of the squares of the light quarks’ charges. The K n and the d n ’s are related to 

the anomalous dimension matrix of the quark and gluon operators [|K] . 

In the absence of QCD, D s ^{x) = ( a em e 2 /2TT ) logQ 2 [1 + (1 — x) 2 ]/x and D 9 ^(x) = 0. 

Substituting this into eq. ED, we recover the singular (in the r —> 0 limit) part of the order 


a em a s result discussed in section II. The relevant part of eq. (|2.3|) is, as promised, 
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It is not important that the x —> 1 limit of eq. (|2.3|) is zero, while that of eq. (|3.3| ) is not. The 


resummed photon spectrum vanishes at x — 1 anyway. Away from x — 1 the resummation of 
the leading logarithms yields finite, non-zero corrections. Although the gluon fragmentation 
into a photon starts only at order a^, its contribution for not too large photon energies is 
comparable to that of the s quark. 

The photon spectrum can be reconstructed from the moments in eq. fl3.2|) by inverse 
Mellin transformation. We use three flavors in the anomalous dimension matrix, as charm 
quarks are excluded from the final state in the experimental analysis. The sensitivity of 
the result to the number of flavors in the /3-function (three or four) is numerically very 
small. The D 9 ^ 1 term, being proportional to the square of the quark charges, would be 
substantially bigger if charm were not excluded. We plot in Fig. 2. the ratio of the resummed 
T 88 contribution to the photon spectrum divided by eq. (|3.3| ). For not too small values of the 
photon energy the resummation suppresses the photon spectrum. However, this suppression 
is not very strong; it only exceeds a factor of two for about x > 0.85. (The spectrum is 
enhanced for x < 0.3.) 

IV. PHENOMENOLOGY AND CONCLUSIONS 

While we think the observations above are interesting even purely from a theoretical 
point of view, we would like to discuss whether there are any experimental implications. 
Although enhanced by log (ml/// 2 ), the r 88 term discussed in this paper cannot become 
arbitrarily large as compared to To- The reason is that in the large mb limit the running 
of a s {rrib) cancels the logarithmic enhancement. For a fixed value of the fi —> 0 limit is 
not allowed. On the one hand, it is cut off by nonperturbative effects; on the other hand, 
the uncalculable parts of the fragmentation functions can be absorbed in /i, as long as /i is 
of order Aqcd- Beyond perturbation theory, the log(l/r) term of the order a em a s result 
should be understood as log (ml/ fi 2 ), where /i is some scale of order Aqcd- Thus, it is also 
obvious that a measurement of the r 88 contribution would not be sensitive to the strange 
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quark mass. 


For our numerical estimates we use C 7 {m b ) ~ —0.31 and C 8 (m b ) ~ —0.15, corresponding 
to m t = 175 GeV and a s (Mz) = 0 . 12 . We also adopt m b = 4.8 GeV and 0.3 < /i< 1 GeV. 
The T 88 contribution to the energy region relevant for CLEO turns out to be negligible. 
The reason is mainly the overall C\f9C 7 suppression. On the other hand, any process that 
contributes to B decays also contributes to B —> X 7 via the fragmentation of any of the 
decay product partons into a photon. We estimate that the tree-level IF-mediated decay 
B —> uud —> X 7 constitutes of order 1% background to one of the CLEO analyses that 
does not reconstruct a kaon. This effect becomes more significant and needs to be taken 
into account more precisely for any future measurements of the 6 —> d 7 decay that aim at 
measuring V t d- 

The T 88 contribution, however, dominates the photon spectrum below about 1.2 GeV. 
Measuring the inclusive B —> X s 7 spectrum at such low photon energies will probably 
remain impossible at CLEO. The situation may be better at asymmetric B factories [|T 6 ] . 
But even if such a measurement is possible in a bin around 1 GeV (say, 0.3 < x < 0.5, that 
corresponds to roughly 0.7 < E 1 < 1.2 GeV), the B —> X s 7 branching fraction into this 
bin is about 6 x ICG'. Taking into account expected experimental efficiencies for such an 
analysis, probably only a handful of events can be detected (assuming 10 8 B' s). The reward 
for this rather hard experimental analysis would be a theoretically clean determination of 
CV 


In conclusion, we pointed out that a part of the next-to-leading order result for the inclu¬ 
sive B —> X s 7 decay rate is proportional to a s log (m b /m%). Summing the leading logarithms 
of the form [a s log(?n^//i 2 )] n to all orders in perturbation theory suppresses (enhances) this 
contribution to the photon spectrum for large (small) values of the photon energy. In the 
large m b limit, formally, the logarithmic factor cancels the a s suppression of the T 88 contri¬ 
bution, as compared to the leading order result To- A measurement of the photon spectrum 
about E 1 ~ 1 GeV would provide a theoretically clean determination of C 8 , the Wilson 
coefficient of the b —*> s g magnetic moment type operator. 
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FIG. 1. Dalitz plot for the contribution of the square of the invariant matrix element a) O 7 , 
and b) Os to the decay b —> s 7 g in the m s —> 0 limit. The shaded area is the available phase 
space, the thick lines denote where collinear singularities occur, and the black box indicates where 
the virtual corrections contribute. 
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Figure 2 
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FIG. 2. The resummed r§8 contribution to the photon spectrum divided by the leading log 
part of the Born term, eq. (3.3), as a function of x = 2E 1 /(1 — r)mf, (solid curve). Dashed curve 
is the gluon decay contribution, eq. (3.2b). 
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